Large contributions to dark matter annihilation from three-body final states 
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The annihilation rate of dark matter particles plays a crucial role in dark matter studies, for it 
determines their relic density and their indirect detection signal. In this paper, we show that this 
annihilation rate can receive large additional contributions from three-body final states consisting of 
a real and a virtual massive particle, such as WW* (—¥ Wff') and ti* (— > tWb). We consider two 
specific examples, from the singlet model and the MSSM, and find that, due to the new three-body 
final state contributions, the prediction for the relic density may decrease by more than a factor two, 
whereas the present dark matter annihilation rate gets enhanced by up to two orders of magnitude. 
Some of the implications of these results are briefly discussed. 



INTRODUCTION 



Dark matter annihilations have two different and 
equally important roles in the phenomenology of dark 
matter models. On the one hand, the annihilation of 
dark matter particles in the early Universe dictates, via 
the Boltzmann equation, the expected relic density and, 
consequently, the viable parameter space of dark matter 
models. This parameter space is usually obtained by re- 
quiring, on top of accelerator and precision bounds, that 
the relic density be consistent with the observed dark 
matter density [1]. Hence, an accurate computation of 
the dark matter annihilation cross section is almost a pre- 
requisite for dark matter studies. In addition, the present 
annihilation rate of dark matter particles determines the 
indirect detection signatures of dark matter. Thus, to 
find the expected fluxes in gamma rays, neutrinos, and 
antimatter generated by dark matter annihilations, we 
must first compute properly the cross section for the an- 
nihilation of dark matter particles. 

Up to now, most studies have only consider dark mat- 
ter annihilations into either two-body final states, or 
three-body final states with a massless gauge boson 
such as qqg or ffj, and have implicitly assumed that 
they dominate the total annihilation cross section. There 
exists, however, another class of three-body final states 
that may have an important impact on dark matter an- 
nihilations. They are final states consisting of a real and 
a virtual massive particle, such as WW* (Wff 1 ) and ti* 
(tWb), and are the main subject of this paper. These 
three-body final states have not been considered in much 
detail in the literature. In fact, neither DarkSUSY Q 
nor micrOMEGAs [4], two sophisticated programs com- 
monly used to study the phenomenology of dark mat- 
ter models, include any of these three-body final states 
in their calculations. To our knowledge, the only works 
where such final states have been considered are [B| and 
@. In Chen and Kamionkowski computed the cross 
section for neutralino annihilation into the three-body 
final states WW* and ti* and showed that, in certain 
regions of the parameter space, they can enhance the 



neutrino signal expected from neutralino annihilations. 
In that work, which has remained largely unnoticed over 
all these years, they did not study the possible effect 
of three-body final states on the neutralino relic abun- 
dance; they simply stated that such effect should be neg- 
ligible. In Q, Hosotani, Ko, and Tanaka considered sta- 
ble biggs bosons as cold dark matter candidates in the 
context of gauge-higgs unification models. They found, 
using a semi-analytic approximation, that the inclusion 
of the three-body final state WW* modifies significantly 
the prediction of the relic density for dark matter masses 
below the W threshold. In this paper, we will update and 
improve these previous results in several ways. We will 
accurately compute the effect of three-body final states 
on the annihilation rate of dark matter particles and on 
their relic abundance in two well-motivated extensions of 
the standard model: the MSSM and the singlet scalar 
model. In both cases, we find that the inclusion of three- 
body final states considerably affects the predicted relic 
density and the indirect detection signal of dark mat- 
ter particles. We thus demonstrate that three-body final 
states consisting of a real and a virtual massive particle 
cannot be generally neglected in dark matter annihila- 
tions. 



THREE-BODY FINAL STATES 

Ordinarily, particle physics processes such as decays 
and annihilations are dominated -provided that they are 
open- by two-body final states. Three-body final states 
are typically suppressed with respect to them by one ad- 
ditional coupling and one more propagator. Exceptions 
to this behaviour are not uncommon though. The decay 
of the higgs boson in the standard model constitute a 
particularly striking example of the importance of three- 
body final states. An intermediate mass higgs boson, 135 
GeV < run < 160 GeV, decays mainly into a three body 
final state consisting of a real W boson and a virtual 
one (W*) that in turn decays into a fermion-antifermion 
pair, h -> WW* -> Wff '% This three-body final 
state accounts for more than 10% of the branching for 



higgs masses above 116 GeV and dominates over bb in a 
wide region below the W + W~ threshold. Thus, to prop- 
erly compute the higgs decay width and its branching 
ratios, two-body decays are not enough; three-body final 
states must necessarily be taken into account. By the 
same token, we will see that, in some cases, to accurately 
compute the dark matter relic density and its indirect 
detection signals it is mandatory to include three-body 
final states. 

Besides illustrating the need to go beyond two-body 
final states, the analogy with the decays of the higgs bo- 
son also allows us to foresee the conditions under which 
three-body final states are expected to be relevant in dark 
matter annihilations. Notice that, in the higgs case, the 
decay h — > W + W~ becomes dominant once it is open, 
and it is below that threshold that the three-body de- 
cay (into WW*) is important. For dark matter anni- 
hilations, accordingly, three-body final states could be 
relevant below the threshold of a standard model parti- 
cle, X, if the annihilation into XX becomes dominant 
for m,DM > mi, where mnu denotes the mass of the 
dark matter particle. Given that, in most dark matter 
models, m DM ~ 100 GeV (WIMPs), X could be a W 
(or Z) boson, a higgs boson, or a top-quark. Hence, 
three-body final states might be relevant for dark mat- 
ter masses below Mw, n^hj or m t . To satisfy the other 
condition observed in the higgs case, we need that the an- 
nihilation into W + W~ , hh or tt become dominant above 
the respective threshold. Thus, whether three-body fi- 
nal states turn out to be relevant or not will depend on 
the specific particle physics model considered and on the 
region of the parameter space examined. 

We will use two examples, singlet scalar dark matter 
and neutralino dark matter in the MSSM, to illustrate 
the relevance of three-body final states in dark matter 
annihilations. In both cases we find regions in the pa- 
rameter space where the inclusion of three-body final 
states modifies in a significant way the dark matter phe- 
nomenology of these models. 




FIG. 1: One of the diagrams with three-body final states 
that enhances the dark matter annihilation rate in the singlet 
model. 




50 60 70 80 

Singlet Mass (GeV) 



FIG. 2: The ratio between the three-body and the two-body 
annihilation rate in the singlet scalar model of dark matter. 



SINGLET SCALAR DARK MATTER 

The singlet scalar model [1, HI is one of the simplest 
extensions of the standard model that may explain the 
dark mattter. It includes one additional scalar field, S, 
that is neutral under the standard model gauge group 
and odd under a new discrete symmetry Zi (S —> —S). 
This symmetry enforces the stability of S and renders it 
a suitable dark matter candidate. This simple model in- 
troduces only two additional parameters to the standard 
model: the singlet mass, ms, and the coupling between 
the singlets and the higgs, A. In addition, the higgs mass 
(m/j), a standard model parameter, also enters into dark 
matter calculations. It has been shown elsewhere, see e.g. 
9], that the singlet scalar model can indeed explain the 
observed dark matter density and that it can be tested 
through direct and indirect detection in future experi- 
ments. 

In this model, singlets annihilate mainly through 
s-channel higgs boson exchange into standard model 
fermions and gauge bosons. Among the two-body fi- 
nal states, the main annihilation channels are bb for 
ms < Mw and W + W~ for ms > Mw- It is a situa- 
tion analogous to that of higgs boson decays mentioned 
previously. Henceforth, it is likely that within the mass 
range 

50 GeV <m s < M w , (1) 

the annihilation cross section is actually dominated by 
the three-body final state WW* (W+W~* + W+*W~). 
Figure [T] shows the Feynmann diagram corresponding to 
the annihilation process into this three-body final state, 
SS -> WW* ->■ Wff, in the singlet model. Notice that 
there are 18 different final states contributing to that 
diagram. How important is this new contribution? 

Figure [2] shows, as a function of the singlet mass, 
the ratio between the annihilation rate, <jv, into the 
three-body final state WW* and that into two-body fi- 
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FIG. 3: The ratio between the relic density obtained including 
the three-body final state and the one predicted for two-body 
final states only. The parameter A was set to 0.03. 
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FIG. 4: The ratio between the relic density obtained including 
the three-body final state and the one predicted for two-body 
final states only. The parameter A was set to 0.03. 



nal states, which is dominated by the bb contribution. 
The figure was obtained for A = 10~ 2 , rr%h = 120 GeV, 
and small v (« 10~ 3 ), typical for dark matter particles 
annihilating in the Galactic halo. To compute the cross 
sections we used the CalcHEP [10| and micrOMEGAs 
Q packages. Even though cru(3-body) and <7t>(2-body) 
are both resonant for ms ~ 60 GeV, their ratio has a 
smooth non-resonant behavior, as observed in the figure. 
Moreover, because A and mu affect both cross sections 
in the same way, the line shown in figure [2] is actually 
independent of these parameters. In other words, it is a 
generic prediction of the singlet model. 

Notice from the figure that already for ms ~ 60 GeV 
the three-body final state is not negligible, accounting 
for about 10% of the cross section. For ms ~ 70 GeV, 
the two-body and the three-body cross sections are of 
the same order. Finally, between mg ~ 70 GeV and 
ms ~ Myy, the singlet annihilation cross section is 
clearly dominated by the three-body final state. In fact, 
close to Myy the cross section into WW* is more than 10 
times larger than the total cross section into two-body fi- 
nal states. Hence, for singlet masses between 60 GeV and 
Mw, the three-body final state cannot be neglected in in- 
direct detection studies, for it modifies considerably the 
total annihilation rate and the branching fractions into 
standard model particles, giving rise to different yields in 
photons, neutrinos, and antimatter. 

Due to the higher dark matter velocity in the early 
Universe and to the thermal averaging of the cross sec- 
tion, the effect of the three-body final state on the sin- 
glet relic density is expected to be smaller. It is signifi- 
cant nonetheless, as illustrated by figures [3][5] Figure [3] 
compares the correct relic density (solid line), obtained 
including two- and three-body final states, with the pre- 
dicted relic density for two-body final states only (dash- 
dotted line), as usually computed in the literature. For 



that figure we set mh = 120 GeV and A = 0.03. To ob- 
tain these results we again used CalcHEP, to compute 
the three-body contribution, and micrOMEGAs, to cal- 
culate the two-body part and to solve the Boltzmann 
equation. Both lines have the same generic behavior. 
The relic density initially increases as the singlet mass 
departs from the higgs resonance (located at ms ~ 60 
GeV), it reaches a maximum, and it then decreases in 
the vicinity of the VF-threshold. The value of the correct 
relic density, however, may differ significantly from that 
predicted for annihilation into two-body final states. At 
ms ~ 75 GeV, for instance, the two-body approximation 
overestimates the cross section by more than a factor 2. 

A similar effect is present also for larger higgs masses, 
as illustrated in figure |4j In this case the relic density 
is a decreasing function of the singlet mass. Again, we 
observe that the two-body approximation overestimates 
the relic density for singlet masses between 50 GeV and 
Mw- A better way to visualize the effect of the three- 
body final state on the relic density is by displaying the 
ratio of the correct relic density (two and three-body) to 
the two-body relic density, as done in figure [5j Notice 
that the ratio is significantly smaller than 1 over a wide 
mass range and that it can decrease even below 0.4, a 
large deviation from the two-body result. From this fig- 
ure we learn that in the singlet model the correct relic 
density could be significantly smaller than the relic den- 
sity obtained for two-body final states. 

A detailed analysis of the resulting new parameter 
space of the singlet model, and of the implications of 
three-body final states for its direct and indirect detec- 
tion prospects will be done in a future work. Next, we 
will see the relevance of three-body final states for neu- 
tralino dark matter. 
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FIG. 5: The ratio between the relic density obtained including 
the three-body final state and the one predicted for two-body 
final states only. The parameter A was set to 0.03. 



NEUTRALINO DARK MATTER 

For neutralino dark matter the most promising region 
where three-body final states may play an important role 
is below the top-quark threshold. That is, for neutralino 
masses in the range 

130 GeV <m x <m t . (2) 

In this mass interval, neutralinos can annihilate into a 
ti* pair followed by the decay of the virtual i into W~b 
(xx ti* tWb) . One of the diagrams that contribute 
to this process in the MSSM is illustrated in figure HI 
For bino-like neutralinos, the dominant two-body anni- 
hilation channel in such mass interval is typically bb. Can 
the three-body final state dominate over bb? 

To ensure that diagrams with three-body final states, 
such as figure El are not suppressed, we will focus on a 
region of the parameter space featuring a relatively light 
stop. Table U shows the supersymmetric parameters, de- 
fined at low energy, that we consider in the following. 
They give rise to models compatible with present bounds 
from accelerator and precision data. The lightest stop, 
which is lighter than all other squarks thanks to the non- 
zero trilinear coupling A t , has a mass of about 260 GeV. 




FIG. 6: One of the diagrams with three-body final states that 
enhances the neutralino annihilation rate in the MSSM. 
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FIG. 7: The ratio between the three-body and the two-body 
annihilation rate in the supersymmetric model defined by Ta- 
ble!] 



The neutralino mass is determined by Mi and is left as 
a free parameter to be varied in the range The re- 
sulting lightest neutralino is always bino-like. 

Figure [7] shows, as a function of the neutralino mass, 
the ratio between the cross section into the three-body 
final state ti* and that one into two-body final states. 
Since the latter is almost constant over the neutralino 
mass range shown, the variation observed in the figure 
is entirely due to the cross section into the three-body 
final state. As before, we took v « 10 -3 and used mi- 
crOMEGAs and CalcHEP to compute the cross sections. 
Notice from the figure that at m x ~ 140 GeV the three- 
body final state already reaches 10% of the total annihi- 
lation cross section. For m x ~ 150 GeV, the three-body 
final state contribute to the cross section as much as the 
two-body states. For larger masses, up to the top thresh- 
old, the neutralino annihilation cross section is entirely 
dominated by the three-body final state. In that region, 
the two-body final state approximation commonly used 
in previous works would fail badly, underestimating the 
neutralino annihilation cross section and misjudging the 
resulting final states. As a consequence, neither the nor- 
malization nor the spectrum of the 7, v, e + and p gener- 
ated in neutralino annihilations would be correct. 

The effect of the three-body final state on the neu- 
tralino relic density is illustrated in figures |8] and HI In the 
calculation, coannihilation effects between the neutralino 
and the lightest stop are automatically taken into ac- 
count. They are not important though, because the mass 
difference between them is significant. Notice that the 
correct neutralino abundance, calculated including the 



Parameter 


Mi, M 3 , m£,fj, 


nig 


A t 


M A 


tan P 


Value 


1 TeV 


500 GeV 


1.4 TeV 


1.5 TeV 


5 



TABLE I: Supersymmetric parameters used. M\ was allowed 
to vary freely. 
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FIG. 8: The ratio between the relic density obtained including 
the three-body final state and the one predicted for two-body 
final states only. 



FIG. 9: The ratio between the relic density obtained including 
the three-body final state and the one predicted for two-body 
final states only. 



three-body final state, could be more than 10% smaller 
than that obtained taking into account only two-body 
final states, as usually done in the literature. 

In view of these results, the claim that DarkSUSY com- 
putes the neutralino relic density, for any MSSM model, 
with an accuracy of 1% ll| no longer holds. To reach 
that accuracy, the three-body final states studied in this 
work must necessarily be included in the calculation. 

In any case, notice that, in contrast with the singlet 
model studied in the previous section, the large effect due 
to three-body final states below the top threshold is not 
a generic feature of the MSSM. It is only within certain 
configurations, for specific sets of parameters, that the 
effect becomes important. It must be also said that there 
certainly exist additional corrections not included in our 
analysis, such as one-loop corrections to the cross sections 
and toponium resonance effects below the top threshold, 
that could also modify the neutralino annihilation cross 
section and the predicted neutralino relic density in the 
MSSM. 

The most significant feature brought about by dark 
matter annihilation into the three-body final state tt* in 
supersymmetric models is that, as we have seen, it may 
give rise to large corrections to the neutralino annihila- 
tion cross section and to the relic density of neutralino 
dark matter. 



CONCLUSIONS 

We have shown that three-body final states consist- 
ing of a real and a virtual massive particle may play 
an important role in dark matter annihilations. Two 
well-motivated scenarios were used to illustrate this ef- 
fect: the singlet scalar and the MSSM. In both cases, we 
demonstrated that these new contributions can signifi- 
cantly modify the prediction for the relic density of dark 



matter and for its indirect detection signatures. In a fu- 
ture work we will study in more detail the implications 
of these results. 
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